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Abstract 


We examine saturation properties of a quark-based picture of nuclear matter. 
Soliton matter consisting of nonlocal confining solitons is used to model nu¬ 
clear matter. Each composite nucleon is described by a non-topological soliton 
as given by the Global Color Model. We apply techniques and concepts from 
the description of crystal lattices. In particular, the Wigner-Seitz approxima¬ 
tion is used to calculate the properties of the soliton lattice at the mean-field 
level. We focus on infinite nuclear matter at around standard nuclear matter 
density with the simplest one-parameter assumption for the gluon propaga¬ 
tor. The saturation density and incompressibility are calculated as functions 
of the single input parameter of the model. 


I. INTRODUCTION 

The Global Color Model (GCM)[]has been developed to model quantum chromodynamics 
(QCD) at energies low on the scale of high-energy particle physics [||]. The successes of the 
GCM in the hadronic sector || include the reproduction of chiral perturbation theory results 
IHHOI], meson form factors |5j and spectra ||, and both, the soliton |7j and the Faddeev |8j 
picture of the nucleon. These successes motivate the application of this particular soliton 
model when a description of nuclear matter in terms of quarks and gluons is attempted. 

Vladimir Naumovich Gribov efficiently used condensed-matter techniques in studies of 
quark confinement [9j. This is but one example of the recent convergence in language and 
methods between nuclear/particle physics and condensed-matter physics. In this paper we 
use the theory of crystal lattices in a description of nuclear matter in terms of quarks and 


1 Not to be confused with the “Generator Coordinate Method” with the same acronym. 


1 



gluons, with the GCM providing the physics of the primitive cell, which corresponds (for 
large values of the lattice constant) to an isolated composite nucleon. 

To model nuclear matter in terms of a many-soliton problem where each individual 
soliton is accountable to the underlying quark and gluon degrees of freedom is an ambitious 
undertaking flT0|-jT^|. At low excitation the kinetic energy of the composite solitons can 
be neglected as a simplification. Due to the short-range repulsion between the solitons, 

ra 


a soliton crystal emerges as the lowest-energy solution at the mean-field level 


This 


crystal, however, is complicated relative to the structures familiar from condensed matter 
physics in that the lattice sites are themselves populated by extended objects (solitons), 
localized to some extent for large lattice constants. As the density increases, we expect a 
mean-field solution which no longer has localized centers of density. 

An early success of condensed-matter physics is the Kronig-Penney model 0], which 
treats a crystal as a one-dimensional periodic set of square wells. This simple description 
allows analytic solutions to be found for the Bloch wave functions and for the energy bands 
usually associated with crystal lattices. In the case of nuclear matter, deeply bound light 
quarks call for a relativistic approach. Glendcnning and Banerjee have solved the relativistic 
Kronig-Penney model and obtained analytic solutions for the eigen-value spectrum which 
retains the band structure 0- 

In more complete geometries analytic solutions of the Dirac equation for the Bloch wave 
functions and eigen-energies are not available, so some kind of approximation is needed. 
Here we apply the Wigner-Seitz approximation JT6| to nuclear matter in the ground state. 
This entails the assumption of a spherically symmetric mean field, which takes into account 
the effect of the surrounding matter (lattice sites) on the primitive cell in a directionally 
averaged manner. This approximation may quite naturally take into account that nuclear 
matter is more like a fluid than a crystal [117 . 


As discussed above, the primitive cell will be given by the GCM, which admits soliton 
solutions with an intrinsically generated, extended qq meson field. In addition, the individual 
GCM solitons are confining, as can be seen by the lack of poles in the quark propagator 
outside the region where the meson field is nonzero [|7|,[l8|. We have recently applied this 
kind of picture to draw attention to the intersection of the ground state band with the next 
unoccupied band at high densities |[T8|.p~9|| . The intersection represents a color insulator-to- 
conductor transition in the model, and thus signals the onset of quark deconfinement. In 
this sense these ideas of condensed matter physics are applied in the spirit of Gribov ||, 
albeit to a model situation in place of full QCD. 

Our goal in the present paper is to look at the saturation properties of this description 
more closely. The paper is organized as follows. In Section 2, we describe the model used 
for the primitive cell. The application of the GCM to nuclear matter is reviewed in Section 
3. In Section 4 we present new results around saturation density, including nuclear-matter 
incompressibility. Section 5 contains the discussion and the summary. 


II. GCM: THE PRIMITIVE CELL 

As quarks and gluons are confined within hadrons, the description of low-energy nuclear 
matter has traditionally been accomplished in terms of effective hadronic degrees of freedom. 
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The effective degrees of freedom should be derived from the QCD Lagrangian. However, until 
an appropriate solution of QCD at low energies and large length-scales becomes available, 
it is necessary to use a model in place of full QCD. Clearly, we attempt to keep as much 
of the essential physics in the model as possible. An approach to systematize the needed 
approximations and to connect QCD to an effective hadronic field theory can be formulated 
in terms of functional integral methods. The strategy is to transform the integration variables 
from quark and gluon fields to hadron fields. 

One particular implementation of the above ideas is in the framework of the GCM, which 
starts with a truncation of QCD [TJ, leading to the Euclidean action^ 


S[qq] = — j d 4 x d 4 y 


q(x)( 7 ' d + m)S(x - y)q{y) + ^-j a (x)D^(x - y)jl{y) 


( 1 ) 


Here j*(x) = q{x)^-'y 1/ q(x) is a local quark current, with Euclidean Dirac matrices and 
Gell-Mann matrices A a . The two-point gluon function, D^, can be considered the phe¬ 
nomenological input point for the model. Using a Feynman-like gauge, D^ u = S /IU D(x — y ), 
the gluon propagator is particularly simple, and provides a single parameter function for the 
GCM. In (D, m is a current quark mass, which will be taken to be zero in the following. 
The GCM has the global symmetries of QCD, but lacks local gauge invariance. Our choice 
of the gluon propagator is dictated by the requirements of simplicity and quark confinement. 
We take a delta function in momentum space, 

9 2 D{q) = 37T 4 a 2 5 {4 \q) . (2) 


The strength parameter a is the only input of the model. This form of the gluon propagator 
has been shown to produce quark confinement and reproduce nucleon properties 


20 


To exhibit nonlocal quark-antiquark structures in the action, a Fierz reordering may be 
performed |2T[ , which transforms the current-current term in (|1|) as 


J d 4 x d 4 y j*(x)D(x - y)j*(y) = ~\j d A x d 4 y J d (x,y)D(x - y)J e (y,x) . 


( 3 ) 


Here, J e (x,y) = q(x)A e q(y) can be looked upon as a quark-antiquark bilocal current with 
quantum numbers specified by 6. The quantity A 6 is a direct product of Dirac, flavor, and 
color matrices, and contains both, color-singlet and color-octet terms. We focus on the 
color-singlet sector in this work, ignoring correlations that correspond to diquark degrees of 
freedom. 

To cast the partition function in terms of Bose fields, auxiliary nonlocal fields, B e (x,y), 
are introduced, and the partition function is multiplied by 


1 


N J VBexp 


d 4 xd 4 y 


B e (x,y)B 0 (y,x ) 
2 g 2 D(x - y) 


( 4 ) 


where N is a normalization constant. After the transformation B e (x, y) —> B e (x, y)+g 2 D(x— 
y)>J e (y,x), the action is bilinear in terms of the quark fields and the Grassman integration 
can be performed in the usual way. This yields the action in terms of bilocal Bose fields. 
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The replacement of the quark fields with Bose fields is, in principle, an exact functional 
change of variables. Observables calculated from the action (|l|) are not affected by the 
variable transformation, but are now expressed in terms of the Bose degrees of freedom, 
provided the entire sum over 6 is kept. This is impossible in practice, and the truncation 
scheme used can be developed into a systematic method of approximation. To retain the 
chiral content of the QCD action, at least two Bose fields need to be kept. 

The classical vacuum configuration B e 0 is identified by 5S/5B 6 = 0. This produces a quark 
self-energy, T,(x — y) = A 6 Bq(x — y ) satisfying a Schwinger-Dyson equation. In momentum 
space 

r \ a 1 \ a 

S(p) = n . P\A(I?) - 1] + Btf) = g*j J^D(P — - 7 , . (5) 


This can be considered an integral equation for the self-energy amplitudes A{jp 2 ) and B(p 2 ). 
Numerical solutions for these amplitudes are now available at varying degrees of sophisti¬ 
cation 0- Our choice for this study is governed by simplicity within the context of the 
requirement of confinement. The delta-function gluon propagator (|2]) simplifies the integral 
equation (|5|) resulting in an algebraic equation. Analytical solutions for A and B can be 
found as |T"S| 


A(v 2 ) -I 2 2 i B(v 2 ) - I ( a2 ~ A P^~ 2 p2 - 

{p) ~ U[i + (i + ¥rt " (p) I o P 2 >- 


A 


( 6 ) 


where the input parameter a controls the coordinate space width of A and B. The lack of 
solutions to the equation p 2 + M 2 (p 2 ) = 0, (where M = B/A is the dynamic quark mass) 
indicates that ([]) produces quark confinement as there is no on-mass-shell point, thus the 
propagation of a quark in the normal vacuum is prohibited JT8[. It is important to note that 
the amplitude B(p 2 ) plays a dual role in the model: it also acts as the distributed vertex 
for coupling the quarks to the qq Goldstone modes ; 7|.|2~i] . 

The fluctuations B e (x, y) = B e (x — y) — B e 0 are identified as the propagating Bose fields. 
If the color-singlet scalar-isoscalar and pseudoscalar-isovector fluctuations are retained, the 
formalism can be adopted to the requirement of chiral invariance by the variable transfor¬ 
mation 


A t 'B e (x,y) = 

J 7T 


(7) 


where r = x — y and R = (x + y)/ 2, are relative and cm-like coordinates, respectively, f n 
is the pion decay constant, and it has been assumed that the on-shell form factor B can 
also be used off-shell. As a further simplification, the 0 = 0 point on the chiral circle can 
be fixed. In this case the radial fluctuations away from the chiral circle coincide with the 
scalar-isoscalar field variable prior to the transformation. In the numerical work that follows 
the single fluctuation field x corresponding to this situation will be used. Letting m —> 0, 
the action (up to a constant) can be written as a sum of fermionic and bosonic terms: 


S\m,x] = ~ Tr \ lnG 1 {V,X) -InG 1 (0,x)] + J d4R [t^x) 2 + U(x 2 )] ■ (8) 
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The chemical potential (/i) dependence of the fermion term in equation ((^) ensures that a 
meson source from the valence quarks will be generated 0. The [/(y 2 ) term is the effective 
meson self-interaction [jlj. For /j, = 0 the inverse quark propagator takes the form 


G 1 ( x .;/ ) = ■ (), A(x - y) + /„ 1 B(x -y)x (-y^) . 


(9) 


and the saddle-point configuration is at y — fn |20 


Since G 1 (x,y) is time independent, stationary eigenstates of the form 'Uj(x) can be 
obtained from a self-consistent Dirac equation, which in coordinate space takes the form 

0 = /d 3 y|(- 74 e i + 7 - V)A(x-y) + S(x-y) + (^y^)}^( y ) • (10) 

The new x — X ~ fn, and as p 4 = g 4 = where €j is the energy eigenvalue, the meson 
vertex B has an energy dependence. It can also be seen that a wave-function renormalization 
appears with the renormalization constant Z. } given by [0 




d 3 p d 3 q u,j (p ) dG P ’ ^ u 3 (q) . 


( 11 ) 


The meson field equation = 0 may be summarized as 


- Vy(z) 


5U 

<$x(z) 


+ Q x ( z) — 0 , 


with the meson source provided by the valence quarks according to 

"x + y 


Q x ( z ) = 517V / d * x d3y u j( x ) b (~^> x ~ y) 5 


— Z 


%(y) 


( 12 ) 


(13) 


Equations ( H(^) and ( H2|) form a system of coupled differential equations for the quark wave 
functions and the meson field, which need to be solved selfconsistently, with the appropriate 
boundary conditions. 


III. INFINITE NUCLEAR MATTER 


The use of soliton matter to model high-density nuclear matter originated in the mid¬ 
eighties Ha- The new feature of the present work is that the GCM soliton represents 
the primitive cell. Some of the above authors also apply the Wigner-Seitz approximation 


H , to which we turn next. 
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FIGURES 


FIG. 1. The top figure shows what the primitive cell looks like for a simplistic two dimensional 
case. The bottom figure shows what the primitive cell looks like in the Wigner-Seitz approximation. 


A. Wigner-Seitz approximation 


As a means of describing nuclear matter, we consider an infinite collection of GCM 
solitons. At the lowest energies the solitons are expected to arrange themselves in a crystal 
lattice |T5| , [T3[| . Accordingly, the single-quark eigen-energies will develop into energy bands. 
For simplicity, we assume a simple cubic crystal (sc). For a periodic lattice, the Dirac wave 
function must be invariant to a lattice translation, so the solutions must have the Bloch 


form 25 


u 


lat i 


r = Um r e 


/m ■ r 


(14) 


where m = (m x , m y , m z ) is the lattice momentum and u m (r) is a Dirac spinor which has 
the periodicity of the lattice. To solve for the Bloch functions we employ the Wigner-Seitz 
approximation [lfj. This amounts to taking the primitive cell as given by the geometry 


of the lattice (a cell bounded by perpendicular bisectors of the lattice vectors as shown in 
Figure |l]) and replacing it by a spherical cell of radius R and solving for m = 0 in (|l4|). 
The full solution «“(r) is then approximated by Wo( r ) e ?m ' r - Changing the density will 
be implemented in our calculation by varying the cell radius R. 

The Wigner-Seitz approximation places specific boundary conditions on the Dirac 
spinors. These conditions express the requirement that the upper component of the wave 
function must be periodic and anti-periodic for the bottom and the top of the band, re¬ 
spectively. We focus attention on the lowest-energy state of the band, for which the above, 
together with the r = 0 boundary conditions, implies 


9'(r)\ r=R = f(r) L = o = 0 , 


(15) 
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where g{r ) and f(r) represent the radial parts of the upper and lower components of the 
Dirac wave function, respectively. In addition, the meson held of equation (|12|), which also 
appears in the source term of the Dirac equation (!1C|) must now be periodic in r, so that 


X(r + 2 R) = x(r) ; x'(r)\ r=R = 0 , 


where y(r) is the radial part of the meson held. 


(16) 


B. GCM on the Lattice 


It is convenient to solve the Dirac equation in momentum space, while the nonlinear 
Klein-Gordon equation is easier to handle in coordinate space. We seek solutions of the 
Dirac equation with the boundary conditions ([15]). These can be incorporated using a three 
dimensional Fourier expansion. Expanding A, B , u, and the meson-held source, we integrate 
over y and use orthonormality to get an equation for the Fourier components of the quark 
wave function 


-eA{k n ) + B{k n ))5 r 
knA{kn)&nm 


k n-^i, k n) ^nm 
(eA(k n ) + B(k n ))5 n 


+ 


4(7T 


J2 k li V o( k n,k r 


771=0 


0 


m =0 


9m 

fm 


= 0, 


where g n = g(k n ), f n = f{k n ), and 


x(k n ~ k m )Pz(cos 9)d(cos 6) . 


To get this hnal form we have written the Dirac wave functions as 


U{kn) = 


g(k„) 

ia ■ k,J(k r 




(17) 


(18) 


(19) 


where k n = — and we used the spherical symmetry of the Wigner-Seitz cell. If m = 

1,2,..., M, then equation ( |T7|) is a 2M-by-2M eigenvalue problem for the energy eigenvalue 
e. The quantity B/A plays the role of a dynamic mass and the scalar part of the self-energy 
B also acts to couple the quarks to the meson held via (|T8|) . The self-energy terms have an e 
dependence which makes this a highly non-linear problem. The V) term (Z = 1, 2) represents 
the Legendre coefficient for the meson held in the presence of the distributed coupling B. 
One needs to solve the Dirac equation (pT|) and the Klein-Gordon equation for the meson 
held flT2|) selfconsistently. 

To solve for a soliton lattice, we hrst pick a starting meson held and search for the lowest 
energy eigenvalue of Eq. ([17]). This means finding the energy e which makes the determinant 
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of equation (|T7|) vanish. We start at e = 0 and work upwards in energy until the determinant 
changes sign. We then use the bisection method to find the root. Care must be taken so 
that the initial steps are sufficiently fine in e not to miss the lowest root. With the root 
in hand, we can solve for the Fourier components of the Dirac wave functions. For this 
we first perform a lower-upper triangular decomposition and use inverse iteration [^]. The 
momentum-space meson-held source term is constructed from the Dirac wave functions and 
we transform the source to coordinate space for use in the nonlinear Klein-Gordon equation 
for the meson held (|T2|). To solve this nonlinear equation, we treat equation (|l^) as a 
functional of y and use Newton’s method. Once the Klein-Gordon equation is solved for 
the new meson held, we start over with the Dirac equation in this modified meson held. We 
iterate until convergence of the quark wave functions is achieved, which takes between three 
to six iterations to reasonable accuracy. 


C. Energy Bands 


As known from condensed-matter physics, each soliton of the lattice contributes one level 
to each energy band [^5|]. In the Wigner-Seitz approximation we need to calculate only the 
energy for the case when the crystal momentum m in equation (|14|) is zero. The Dirac-Bloch 


) is zero. 

wave function, in the Wigner-Seitz approximation, for an arbitrary crystal momentum is 


«m ( r ) = “o(r) e 


/m • r 


( 20 ) 


and we use fl20|) to calculate the expectation value for the square of the Dirac Hamiltonian 
to estimate the lattice momentum dependence of the energy levels in the band as 


e m = e 


bot 


+ m 


( 21 ) 


where e^ot is the energy of the bottom band. 

To obtain the possible values of m we refer back to the underlying lattice structure. For 
a simple cubic crystal of N solitons and sides of length L = 2RN , the allowed values of the 
component of the lattice momentum in the direction of any of the three axes are 


2tt 47t Nit 



( 22 ) 


with the top of the energy band corresponding to m 
band we obtain 


Ntt 

~L 


7T 

2R' 


Thus for the top energy 


£top 



(23) 


We have performed the same estimate assuming body-centered and face-centered cubic 
lattices. This variation in assumed crystal structure introduces a less than 10 % uncertainty 
in our results. 

The lowest band (l — 0, l' — 1 and j = |) is labeled lsi/ 2 . The next lowest (at 
around standard nuclear density) band has nonzero orbital angular momentum in the large 
















component of the Dirac wave function (/ = 1, l' = 2, and j = §) and is labeled lp 3 / 2 . The 
next band is again an s-state, corresponding to a radial excitation. For very low density 
(R —> cx)) the energy bands shrink to single levels and in the limit reproduce the energies of 
a single soliton. As the density increases, the bands spread out and approach each other. The 
intersection of the energy bands as the density is increased beyond twice standard nuclear 
density is discussed elsewhere [|HJ. However, properties of the model at around saturation 


density have not been examined in detail in Ref. |H|. This is the goal of the present study. 


IV. RESULTS 


Here we focus on the properties of the lowest energy band (ground state band) at around 
standard nuclear density, p « po = 0.17 frn -3 . With one soliton in each Wigner-Seitz cell, 
the energy of each primitive cell is given by the soliton energy. As mentioned above, the 
kinetic energy of the center-of-mass of the soliton is neglected in the present approximation. 
The total soliton energy has a quark contribution which contains the quark kinetic and 
potential energies. Assuming three-fold degeneracy of the energy levels (three colors), the 
quark energy thus equals three times the lowest quark eigen-energy (e) in the self-consistent 
meson held. In addition, the kinetic and potential energy associated with the meson held 
contribute to the soliton energy. The meson held has a kinetic energy given by 


KE X 




The potential energy of the meson held is given by 


PE X = 


fd 3 xU ( X ) 


(24) 


(25) 


where U(x) is the meson self-interact ion from eq. (Sj), expressed as 
tuation held x — X ~ fir- U(x) has the usual “Mexican hat” shape 
takes the form 


a function of the huc- 
M - The total energy 


E = 3e + KE X + PE X 


(26) 


This quantity is plotted against the density in terms of standard nuclear density in Figs. 
2-4 for different values of the parameter a. We selected a range for the values of a earlier 
0 • It was argued that while reproducing the pion decay constant f n closely is particularly 


important for studies at the hadron level, for the present work, which is concerned with the 
density of nuclear matter, the root mean square (rms) radius of the proton constitutes the 
most sensible benchmark. It is not realistic to expect that this simplified one-parameter 
model, which does not explicitly take the pion held into account, should simultaneously fit 
both of these quantities. Since the explicit pion cloud is neglected in the model and it is 
normally assumed that the pions will increase the value of the rms proton radius by about 
20-30 %, we need to choose a to give an rms proton radius ~ .6 or .7 fm. We are able 
to achieve this for a ~ 1.35 ± .1 GeV. Thus, we display calculated results for the values 
a = 1.25,1.35 and 1.45 GeV in the following. 
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In Fig. 2 we show the quark contribution and the total soliton energy for a = 1.25 
GeV. On one of the calculated points we indicate a typical uncertainty we associate with 
our computation. The main source of this uncertainty, which is not more than a couple 
percent, is the freedom in prescribed tolerances at different stages of the calculation. We see 
that as the density increases from zero, the ground-state energy develops a minimum. The 
low-density attraction between the solitons is a consequence of the boundary conditions on 
the quark wave functions 


FIG. 2. The total energy E and the quark contribution to the total energy in MeV, as a function 
of p/po for a = 1.25 GeV. 

In particular, the upper component of the quark wave function is forced to have less 
curvature than in the case of a single soliton, leading to a lower value of the quark kinetic 
energy. At higher densities, where the solitons and the quark wave functions begin to overlap, 
the resulting repulsion overcomes the attraction and the ground-state quark energy starts 
to increase. For this value of a there is a minimum in the total energy at p m in/Po = 0.66. 
To calculate the minimum we performed a second-order fit to the calculated points around 
the minimum, up to p/po = 1.5. The value of the total energy at the minimum is 1844 
MeV. However, our scheme does not take into account the center-of-mass kinetic energy 
associated with the quarks. An estimate of the soliton mass (on the lattice) can be obtained 
by subtracting this motion; M s = [E 2 — (P 2 )] 1 / 2 , where (P 2 ) = 3 (p 2 ), with (p 2 ) being the 
expectation value of the square of the quark momentum. A numerical estimate along these 
lines leads to a soliton mass ~ 15 % below the soliton energy, yielding « 1560 MeV for the 
soliton mass for a = 1.25 GeV. Scaling back to a = 1.05 GeV (used in the earlier single- 
soliton work) would decrease this by another 20 % to ~ 1,200 GeV, not too far from the 
average mass of the nucleon-delta system. Explicitly including the pion will lower this value 
even further |20| . 

Figs. 3 and 4 contain the same information for a = 1.35 and 1.45 GeV, respectively. In 
Fig. 3 we observe a minimum of the total energy at pmin/po = 0.73. Fig. 4 yields a minimum 
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at pmin/Po = 0.98. We have included all calculated points below p/p 0 = 1.5 in all these fits. 
(We carried out calculations at higher densities for each value of a in connection with our 
studies of the QGP transition |L9|. However, in the present work we focus on the second- 
order behavior of the energy versus density curve near its equilibrium value (minimum). We 
estimate an uncertainty in the position of the minima on the order of 15%, as a consequence 
of combining uncertainties in crystal structure, tolerances and fits. 


FIG. 3. The total energy E and the quark contribution to the total energy in MeV, as a function 
of p/po for a = 1.35 GeV. 

Finally we calculate the incompressibility defined as 

dE 

b (p Pmin ) 9p — 

op 

around the observed minima. The resulting values are summarized in Table 1, together with 
the pion decay constant, the root mean square proton (single soliton) radius and the value 
of the density where the minimum occurs for each value of a. As a increases, the root mean 
square radius of a single soliton is moving away from the experimental value of the rms radius 
of the proton, (r 2 )z ~ 0.83 fm. However, as discussed above, the effect of the pion cloud is 
expected to increase the root mean square radius by about a third. The incompressibility 
appears to be in the generally accepted range for the larger values of a. The ~ 20 % 
uncertainty in this quantity is again due to a combination of crystal-structure, tolerance 
and fit-induced uncertainties. As mentioned above, the experimental value of the pion 
decay constant, f n = 93 MeV is best reproduced in our simple single-parameter description 
for a = 1.05 GeV. The value of f n is important for calculations in the hadronic sector. It 
is anticipated that with a more complicated gluon propagator with several parameters, the 
quantities f n , p m in and K can all have reasonable values simultaneously. 


(27) 
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FIG. 4. The total energy E and the quark contribution to the total energy in MeV, as a function 
of p/po for a = 1.45 GeV. 


V. DISCUSSION 


We have calculated saturation properties of nuclear matter with a nonlocal confining 
soliton model, the Global Color Model (GCM). The encouraging successes of the GCM on 
the hadronic level triggered the use of the model for the description of nuclear matter in 
terms of the underlying quark and gluon degrees of freedom. Since there is a large number 
of degrees of freedom to deal with in each individual soliton, the soliton matter picture 
of nuclear matter is necessarily restricted to the mean field level and the Wigner-Seitz 
approximation is utilized. This model provided robust conclusions on the QGP transition 


at high density |T3JT^|. We therefore felt it necessary to investigate the properties of the 
model at around standard nuclear matter density. 

A reasonable reproduction of nuclear matter saturation properties has been found. Some 
caution is necessary though when interpreting the physical significance of this agreement, 
since the kinetic energy of the soliton as a whole has been neglected in our calculation, 
while the Fermi kinetic energy of the nucleons plays an important role in nuclear saturation. 
Another reason while the agreement does not signal a deep physical understanding of nuclear 
matter at the quark level is the fact that we used an extremely simplified, “bare-bones” one- 
parameter gluon propagator. Furthermore, explicit pions were not included in the model. 
We find it remarkable that such a simple effective description nevertheless gives values of the 
equilibrium density and of the incompressibility in the right ball park. We wish to inform 
the community about the existence of such a powerful effective description. 
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TABLES 


TABLE I. Variations of the pion decay constant, the root mean square radius of a single soliton, 
the density where the total energy E is minimum, and the incompressibility with the single model 
parameter a. 


a (GeV) 

1.25 

1.35 

1.45 

U (MeV) 

111 

120 

129 

(r 2 ) 1 / 2 (fm) 

.71 

.67 

.64 

Pmin/ PO 

0.66 ±0.10 

0.73 ±0.11 

0.98 ±0.15 

K (MeV) 

137 ± 30 

220 ± 40 

309 ± 50 


It is of interest to calculate in-medium properties as a function of the density with the 
model. We made the first steps in this direction elsewhere [27] by calculating the axial-vector 
coupling constant and correlation function in medium. It should be kept in mind that more 
realistic gluon propagators are now available in the literature J^j, and should improve the 
model and accommodate a simultaneous fit of the pion decay constant and the root mean 
square radius of the proton. A more complete model requires the treatment of explicit pion 
degrees of freedom for inclusion into the calculation of the root mean square proton radius. 
Developments along these lines will also facilitate the exposure of the chiral content of the 
model and direct comparisons to QCD-based calculations. With such improvements the 
model can become a practical tool in the description of high-density strongly interacting 
matter. 
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